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1 Introduction 

Let <&(x;y) be a meromorphic function on C m x C™ in the variables x = (x\ , . . . , x m ) and 
y = (yi, . . . ,y n ), and consider two operators A x , B y which act on meromorphic functions in x 
and y, respectively. We say that $>(x;y) is a kernel function for the pair (A x ,B y ) if it satisfies 
a functional equation of the form 

A x ®{x;y) =B y ®(x;y). 

In the theory of Jack and Macdonald polynomials [13], certain explicit kernel functions play 
crucial roles in eigenfunction expansions and integral representations. Recently, kernel functions 
in this sense have been studied systematically by Langmann |12| [T3] in the analysis of eigenfunc- 
tions for the elliptic quantum integrable systems of Calogero-Moser type, and by Ruijsenaars 
[2 H 1221 [23] for the relativistic elliptic quantum integrable systems of Ruijsenaars-Schneider type. 

In this paper we investigate two kinds of kernel functions, of Cauchy type and of dual Cauchy 
type, which intertwine pairs of Ruijsenaars difference operators. In the cases of elliptic difference 
operators, kernel functions of Cauchy type for the (A n -i, A n -i) and (BC n , BC n ) cases were 
found by Ruijsenaars [22\ 123]. Extending his result, we present kernel functions of Cauchy type, 
as well as those of dual Cauchy type, for the (BC m , BC n ) cases with arbitrary m, n (under certain 
balancing conditions on the parameters in the elliptic cases). For the trigonometric difference 
operators of type A, kernel functions both of Cauchy and dual Cauchy types were already 
discussed by Macdonald [H] . Kernel functions of dual Cauchy type for the trigonometric BC n 
cases are due to Mimachi [15] . In this paper we develop a unified approach to kernel functions 
for Ruijsenaars operators of type A and of type BC, with rational, trigonometric and elliptic 
coefficients, so as to cover all these known examples in the difference cases. We expect that our 
framework could be effectively applied to the study of eigenfunctions for difference operators of 
Ruijsenaars type. 

As such an application of kernel functions in the trigonometric BC cases, we derive new 
explicit formulas for Koornwinder polynomials attached to single columns and single rows. This 
provides with a direct construction of those special cases of the binomial expansion formula for 
the Koornwinder polynomials as studied by Okounkov [18] and Rains [19]. We also remark that, 
regarding explicit formulas for Macdonald polynomials attached to single rows of type B, C, D, 
some conjectures have been proposed by Lassalle [TT]. The relationship between his conjectures 
and our kernel functions will be discussed in a separate paper. 

Our main results on the kernel functions for Ruijsenaars difference operators of type A and 
of type BC will be formulated in Section [2 In Section [3] we give a unified proof for them on the 
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basis of two key identities. After giving remarks in Section U] on the passage to the g-difference 
operators of Macdonald and Koornwinder, as an application of our approach we present in 
Section [5] explicit formulas for Koornwinder polynomials attached to single columns and single 
rows. We also include three sections in Appendix. We will give some remarks in Appendix [A] 
on higher order difference operators, and in Appendix [B] make an explicit comparison of our 
kernel functions in the elliptic cases with those constructed by Ruijsenaars [23J. In Appendix ICl 
we will give a proof of the fact that certain Laurent polynomials, which appear in our explicit 
formulas for Koornwinder polynomials attached to single columns and single rows, are special 
cases of the BC m interpolation polynomials of Okounkov [18]. 

2 Kernel functions for Ruijsenaars operators 
2.1 Variations of the gamma function 

In order to specify the class of Ruijsenaars operators which we shall discuss below, by the 
symbol [u] we denote a nonzero entire function in one variable u, satisfying the following Riemann 
relation: 

[x + u] [x — u] [y + v] [y — v] — [x + v] [x — v] [y + u] [y — u] 

= [x + y][x — y][u + v][u — v] (2.1) 

for any x, y,u,v £ C. Under this condition, it is known that [u] must be an odd function 
([— u] = —[u], and hence [0] = 0), and the set Q of all zeros of [u] forms an additive subgroup 
of C. Furthermore such functions are classified into the three categories, rational, trigonometric 
and elliptic, according to the rank of the additive subgroup J] C C of all zeros of [u]. In fact, 
up to multiplication by nonzero constants, [u] coincides with one of the following three types of 
functions: 

(0) rational case: e (cm 2 ) u (17 = 0), 

(1) trigonometric case: e (aw 2 ) su\{ttu / u)\) (0 = 

(2) elliptic case: e (an 2 ) cr(u; Q) (ft = Za>i © Za^), 

where a S C and e(u) = exp(2ir\/—lu). Also, o~(x; denotes the Weierstrass sigma function 

(l_^) e ^+^ 

associated with the period lattice Q = Za>i © Zu;2, generated by uj\, uj2 which are linearly 
independent over R. 

We start with some remarks on gamma functions associated with the function [u]. Fixing 
a nonzero scaling constant 5 G C, suppose that a nonzero meromorphic function G(u\S) on C 
satisfies the difference equation 

G(u + S\5) = [u]G(u\5) (ueC). (2.2) 

Such a function G(u\8), determined up to multiplication by 5-periodic functions, will be called 
a gamma function for [uj. We give typical examples of gamma functions in this sense for the 
rational, trigonometric and elliptic cases. 

(0) Rational case: For any 5 £ C (5 ^ 0), the meromorphic function 

G(u\S) = 5 u l s T{u/5), 

defined with the Euler gamma function T(u), is a gamma function for [u] = u. 
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(1) Trigonometric case: We set z=e(u/u)\) and q=e(5/ui), and suppose that Im(5/(Ji)>0 
so that l^l < 1. We now consider the function 



i i 



[u] = Z2 — z 2 = 2y — 1 sm(-Ku/uJi). 
For this [u] two meromorphic functions 

(f )) 



G.(u\5) = { , G + (u\S)=e(^- i (% 5 ))( q /z;qU (2.3) 



OO 

where (z; q)^ = Y\ (1 — q l z), satisfy the difference equations 

G-(u + 5\5) = -[u]G-(u\5), G+(u + S\S) = [u]G+(u\S), 

respectively. Namely, for e = ±, G e (u\5) is a gamma function for e[u]. (Note that the quadratic 
function Q) = \u{u — 1) satisfies ("^ 1 ) = (2) + u -) Another set of gamma functions for e[u] of 
this case is given by 

G+(u\S) = ^ , \ Jo ° , G-(u\S) 



by means of the infinite products with base ga. 

(2) Elliptic case: Let p, q be nonzero complex numbers with \p\ < 1, \q\ < 1. Then the 

Ruijsenaars elliptic gamma function 

T(z;p,q) = y iP f )ao , (z-^qU = f[ (l-pV*) 
{z,p,q)oo (J=0 

satisfies the g-difference equation 

F(qz;p,q) = 6(z;p)T(z;p,q), 0(z;p) = (2; J>)oo(p/2;p)oo- 

Note also F(pq/z;p,q) = Y{z;p,q)~ l . We set p = e^/wi), q = e{5/ui\) and z = e{u/ui), and 
suppose that Im^/^i) > 0, Im(5/wi) > so that \p\ < 1, \q\ < 1. Instead of 0{z;p) above, we 
consider the function 

1 



[u] = -z 20(z;p) = 2V-lsm(iru/u;i)(pz;p) 00 (p/z;p) 00 

which is a constant multiple of the odd Jacobi theta function with modulus so that [u] 

should satisfy the Riemann relation. Then the two meromorphic functions 



OO 

00 



G-(u\5) = e( -^))T(z;p,q) = e( - ^))^^ 

G + (u\6) = e(£Cff))lW,g) = e(&(%> )) ■ ( 2A ) 

satisfy the difference equations 

G_(u + S\S) = -[u]G-(u\S), G+(u + S\S) = [u]G+(u\S), 
respectively. Another set of gamma functions for e[u] is given by 

G e (u\S) = 0(ez^;q^)T(z;p, q) = T (ep^ z^ ; p^ , q^)F ( - ez^;p^,q^) (e = ±). 

In the limit as p — > 0, these examples recover the previous ones in the trigonometric case. 

We remark that, if G(u\5) is a gamma function for [u], then G(5 — "u|<5) _1 is a gamma function 
for —[u]. Also, when we transform [u] to [u]' = ce(au 2 )[u], a gamma function for [u]' is obtained 
for instance as G'(u\5) = c u / s e(ap(u))G(u\5), where p(u) = -^u 3 — ^u 2 + |u. 
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2.2 Kernel functions of type A 

The elliptic difference operator which we will discuss below was introduced by Ruijsenaars |20| 
together with the commuting family of higher order difference operators. In order to deal with 
rational, trigonometric and elliptic cases in a unified manner, we formulate our results for this 
class of (first order) difference operators in terms of an arbitrary function [u] satisfying the 
Riemann relation. (As to the commuting family of higher order difference operators, we will 
give some remarks later in Appendix lAl) 

Fix a nonzero entire function [it] satisfying the Riemann relation (|2.1|) . For type A, we 
consider the difference operator 

m r 1 

»,i.r:M7;;: ; . A i (x;n) = J] & + * J 

. 1 , \Xi—Xj\ 

i=l l<j<rn\j^Fi 

in in variables x = (x\, . . . , x m ), where 5, k £ C are complex parameters ^ f2, and stands for 
the 5-shift operator 

T xJ( x ^ ■ ■ ■ i x m) = f(xi, ...,Xi + 5,... ,x m ) (i = 1, ...,m). 

Note that this operator remains invariant if one replaces the function [it] with its multiple by 
any nonzero constant. 

By taking a gamma function G(u\5) for (any constant multiple of) [it] as in (|2.2p . we define 
a function &a(,%', y\S, k) by 

$ a i;#,k = I 1—^7 — - — - . x . (2.5) 

with an extra parameter v. We also consider the function 

m n 

M«ij)=nii[ j! j-tt+i']' M 

j=lZ=l 

These two functions &a( x ', y\S, k) and ^^(^5 ?/) are kernel functions o/ Cauchy type and o/ tiixal 
Cauchy type for this case, respectively. 

Theorem 2.1. Let [it] be any nonzero entire function satisfying the Riemann relation. 

(1) If m = n, then the function &a( x 'i y\S, k) defined as <\2.5\) satisfies the functional equation 

D^ K ^ A (x; y\S, k) = D^ A (x; y\5, «). 

(2) Under the balancing condition mn + n5 = 0, the function ^a{x;v) defined as (|2.6p satisfies 
the functional equation 

[K]D^ A (x;y) + [S\D^ A {x;y) = 0. 

Statement (1) of Theorem 12.11 is due to Ruijsenaars [22\ 123] . (See Appendix IB. II for an 
explicit comparison between <& A {x\y) and Ruijsenaars' kernel function of [23].) In the scope of 
the present paper, the balancing conditions (m = n in (1), and ran + n5 = in (2)) seem to be 
essential in the elliptic cases. In the context of elliptic differential operators of Calogero-Moser 
type, however, Langmann [13] has found a natural generalization of the kernel identities of 
Cauchy type, which include the differentiation with respect to the elliptic modulus, to arbitrary 
pair (m, n). It would be a intriguing problem to find a generalization of this direction for elliptic 
difference operators of Ruijsenaars type. 

In trigonometric and rational cases, these functions $ A (x;y\5, k) and ^A(x',y) satisfy more 
general functional equations without balancing conditions. 
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Theorem 2.2. Suppose that [u] is a constant multiple of s'm(iru/uJi) or u. 

(1) For arbitrary m and n, the function & A (x; y\S, k) satisfies the functional equation 

D^^ A (x;y\S,K) - D^ a (x;v\S,k) = [{m 7 <Ms; y\6, k). 

[K\ 

(2) The function *$> A (x;y) satisfies the functional equation 

[K]D<i^ A {x-y) + [5}D^ 5 ^ A (x;y) = [mn + nS}V A (x; y). 

These results for the trigonometric (and rational) cases are essentially contained in the dis- 
cussion of Macdonald |14] , 

A unified proof of Theorems 12.11 and 12.21 will be given in Section [3j We will also explain in 
Section H] how Theorem 12.21 is related with the theory of Macdonald polynomials. 

2.3 Kernel functions of type BC 

The (first-order) elliptic difference operator of type BC was first proposed by van Diejen [2]. It 
is also known by Komori-Hikami [9] that it admits a commuting family of higher order difference 
operators. In the following we use the expression of the first-order difference operator due to [9], 
with modification in terms of [it]. (In Appendix iBj we will give some remarks on the comparison 
of our difference operator with other expressions in the literature.) 
For type BC, we consider difference operators of the form 

m m 

E m«) =J2At(x;»\6,K)Tl +^2Ar( x ^\5,K)T-f + A°( X ^\5,K) (2.7) 

i=l t=l 

including 2p parameters /i = (fi±, . . . , fj,2 P ) besides (5,k), where p = 1, 2 or 4 according as 
rankfi = 0, 1, or 2. In the trigonometric and rational cases of type BC, we assume that the 
function [it] does not contain exponential factors. Namely, we assume that [it] is a constant 
multiple of one of the functions 

u (p = l), sm(TTu/ui) (p = 2), e(ait 2 ) a(u;Q) (p = 4). (2.8) 

In each case we define u±, . . . , iv p S as 

(0) rational case: O = 0, u± = 0, 

(1) trigonometric case: f2 = 7Luj\, 102 = 0, 

(2) elliptic case: £1 = Zu>i © Zlo2, ^3 = — oj\ — u>2, loa = 0. 

Then the quasi-periodicity of [it] is described as 

[it + u r ] = e r e(r) r (u + ~a; r )) [u] (r = 1, . . . , p) 

for some rj r G C and e r = ±1. In the trigonometric and rational cases (without exponential 
factors), one can simply take r] r = (r = 1, . . . , p). Note also that [it] admits the duplication 
formula of the form 



[2 „j=2 H n4^i, p^-n ^te^ . (2.9) 
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(In the trigonometric and rational cases, these formulas fail for [u] containing nontrivial expo- 
nential factors e(au 2 ).) In relation to the parameters \i = (p,\, . . . ,p>2p), we introduce 



8=1 S=l 



U3 K . 



This parameter is related to quasi-periodicity of the coefficients of the Ruijsenaars operator, 

p 

and plays a crucial role in various places of our argument. Note that the last term uj s is 

8=1 

nontrivial only in the trigonometric case: it is uj\ when p = 2, and when p = 1, 4. 

With these data, we define the coefficients of the Ruijsenaars operator E^ S ' K ^ of type BC 
as follows: 

2p 

nfri + M.] r r . + T ., K i 

At(x;p J \6,K) = - H 



[si - + ±(<J - uj s )} i<j<m;j& ^ ± Xj 



s=l 



A i (x;ii\6,K)=Af(-x',ii\5,K) (i = l,...,m), (2.10) 

with an abbreviated notation [x ± y] = [x + — y] of products, and 

p 

A°(x; fi\5, k) = ^^A°(x; p,\S, k), 

r=l 

e{-(mK + ±cMW)Th.)n$(u>r-S)+t*.] ™. , 1( , 

m n tiK - «.)] n [|k - ^ + * - *)]i=i 1 dj ± Xjl 

Sj^r s=l 

This operator E^ 5,K ^ is one of the expressions for the first-order Ruijsenaars operator of type 
BC m due to Komori-Hikami [91 (4.21)]. We remark that this difference operator Ex has 
symmetry 

E M-S,-k) = E (-h\8,k) (2.12) 

with respect to the sign change of parameters. (This property can be verified directly by using 

p 

the quasi-periodicity of [u] and the fact that Yl e s = 1, — 1, — 1 according as p = 1,2,4.) Note 

8=1 

also that it remains invariant if one replace the function [u] with its multiple by a nonzero 
constant. By the duplication formula (|2.9p . one can also rewrite the operator E X ^ S in the 
form 

2p 

p-l m I! [ x i + Ms] r _i_ . I 

iTTrL , i 2 p(mI<5a) _ V_f=l TT <Ei Xj + KJ r& 

-MI- 2 u> s \e x xi + §] 11 



s=l i=l L 111 1 1 l<j<m;j^i 

m Tl i~ x i +Vs] r , , i 

, \ ^ 8 = 1 TT [— x i ± x j T ,5 

+ ^[-2x i ][-2x i + fl 1 ^ll^ [-x^x,] - 
where i<T r = e( - (w r + (m + 1)k - 5 + \M 6 ^)rj r ). 



(2.13) 
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By using any gamma function G(u\5) for (a constant multiple of) [u], we define a function 
&bc(x; y\5, k) of Cauchy type, either by 

fJit^Gixj ±yi + %(8 + K)\8) 

or by 

m n 

$ BC (x;y\5,K) = JJJJ JJ + e 2 y z + ±(5 - k)|5). (2.15) 

j=ll=l ei,62=± 

Note that the function (|2. 15j) differs from (|2.14j) . only by a multiplicative factor which is 5- 
periodic in all the variables Xj (j = l,...,m) and yi (I = l,...,n). We also introduce the 
function 

m n 

*Bc(x;y) = l[l[[x j ±y l } (2.16) 

j=ll=l 

of dual Cauchy type. 

Theorem 2.3. Suppose that [u] is a constant multiple of any function in (j2.8j) . 

(1) Under the balancing condition 

2(m - n)« + c^' 5 ' k) = 0, 
£/te function ^Bc( x ',y\S, n) defined as (|2.14p or (|2.15j) satisfies the functional equation 

EM s ^^ BC (x;y\8,K) = E^$ B c(x;y\5,K), 
where the parameters v = (i/i, . . . , V2 P ) for the y variables are defined by 

v r = -(<$ + k) - fi r (r = 1, . . . , 2p). 

(2) Under the balancing condition 

2mn + 2n5 + c M5 ' k) = 0, 
the function ^Bc( x ',y) defined as (|2.16p satisfies the functional equation 

[K]^IW« BO (i;y) + [<S]^ m *bc(*;i/) = o. 

Statement (1) for the cases m = n is due to Ruijsenaars \21\ I23j . An explicit comparison 
will be made in Appendix IB. 21 between our $Bc(x;y\5, k) and Ruijsenaars' kernel function of 
type BC in [23]. 

In the case of variables, the Ruijsenaars operator E^ S,K ^ with m = reduces to the 
multiplication operator by the constant 



q(h\6,k) _ 



e (_i c MW^) n[iK-<5)+M-] 



L2 1 



sj^r s=l 
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We remark that Theorem 1 2, 3 \ (2) for n = implies that, when 2mn + c^ s ' K ' = 0, the constant 
function 1 is an eigenfunction for E x K ^'- 

(Statement (1) for n = gives the same formula since C^ S ' K ^ = — [S]C^ K ' S ^ /[k].) 

In the trigonometric and rational cases, these functions <&bc{ x \ y\^i k ) an d ^Bc( x ]y) satisfy 
the following functional equations without balancing conditions. 

Theorem 2.4. Suppose that [u] is a constant multiple of s'm(iru/uJi) (p = 2) or u (p = 1). 

(1) The function $>bc( x 'i y\S, k) satisfies the functional equation 

[n]E^ s ' K ^Bc(x;y\5,K) - [K]E^ s ' K ^Bc(x;y\S,K) 
= [2(m - n)K + c^ s ^}$ BC (x; y\8, k), 

where v r = ^(5 + k) — p r (r = 1, . . . , 2p). 

(2) The function ^Bc( x ',y) satisfies the functional equation 

[ K ]E^ 8 > K ^Bc(x;y) + [8}E^ K > 5 ^Bc(x;y) = [2mn + 2nd + C W s >«)]y BC ( X ; y). 

As a special case n = of this theorem, we see that the constant function 1 is a eigenfunction 
of E X ^ S,K ^ for arbitrary values of the parameters. Theorems 12.31 and 12.41 will be proved in 
Section 

In the trigonometric and rational BC cases, it is convenient to introduce another difference 
operator 

m m 

D W, K ) = ^2A+(x;p\5,K){Tl - 1) + J2A7(x;p,\6, K )(T-* - l) 

i=l i=l 

with the same coefficients A\{x\ p\5, k) (i = 1, . . . , m; e = ±) as those of the Ruijsenaars opera- 
tor E X ^ S ■ In the trigonometric case, this operator D X ^ S,K ^ is a constant multiple of Koorn- 
winder's (/-difference operator expressed in terms of additive variables. By using the relation 

one can rewrite the functional equations in Theorem 12.41 into those for D X ^ S ' K \ 
Theorem 2.5. Suppose that [u] = 2\J — 1 sin(7m/a>i), p = 2 [resp. [u] = u, p = 1). 

(1) The function $>bc( x 'i y\S, k) satisfies the functional equation 

[ K ]DMW$Bc( x ;y\6,K) - [k]d^^bc( x ;v\6,k) 

= [m/t][- nK][(m — u)k + c^ S ' k ^]^bc{ x ', y\S, k) (resp. = 0), 

where v r = | (5 + k) — p, r (r = 1, . . . ,2p). 

(2) For arbitrary m, n, the function ^>Bc{ x ]y) satisfies the functional equation 

[K]D^ 5 ^ BC ( x ;y) + [5]D^ 5 ^ BC ( x ;y) 

= [mK][n5][mK + n5 + c^ s ^}^ Bc{ x ]y) [resp. = 0). 

Statement (2) of Theorem 12 . 5 1 recovers a main result of Mimachi [15]. A proof of Theorem 12. 51 
will be given in Section 13.41 Also, we will explain in Section H] how Theorem 12.51 is related to 
the theory of Koornwinder polynomials. 
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3 Derivation of kernel functions 

3.1 Key identities 



We start with a functional identity which decomposes a product of functions expressed by [u] 
into partial fractions. 

Proposition 3.1. Let [u] be any nonzero entire function satisfying the Riemann relation. Then, 
for variables z, (x\, . . . , xn), and parameters (ci, . . . , cat), we have 

j=l 1 jJ i=l 1 li l<j<N;j& 1 1 jl 

where c = c\ + • • • + c^r. 

This identity can be proved by the induction on the number of factors by using the Riemann 
relation for [u]. Note that, by setting yj = xj — Cj, this identity can also be rewritten as 

V . l z ~ Vj\ _ V^r„, „.A z ~ x i + °] TT ~ 



11 Z-X, \Z — Xj\ xx 

j=l L n i=l L 11 Kj<N: 



<j<N;j^ [Xl X ^ 
N 

where c = Yl ( x j ~ Uj)- 

3=1 

From this proposition, we obtain the following lemma, which provides key identities for our 
proof of kernel relations. 

Lemma 3.1. Consider N variables x±,... ,xn and N complex parameters c\, . . . , cjv- 

(1) Let [u] be any nonzero function satisfying the Riemann relation. Then under the balancing 

N 

condition c j = 0, we have an identity 

3=1 

N 

En II i =0 - ( 3 - 2 ) 

^-^ ±A - \Xi — Xj\ 

*=1 l<3<N;jy£i 1 n 

as a meromorphic function in (x%, . . . ,xn). 

(2) Suppose that [u] is a constant multiple of sin(-7ru/u;i) or u. Then for any c±, . . . , c/y £ C, 
we have an identity 

N 



N 

EM II [x \~-ti 

i=l Kj<N;M [l jJ 



l<j<N;tf 

as a meromorphic function in (x\, . . . , xjy 



i=l . 



(3.3) 



When [z] = sin(7rz/6Ji), we have [z + a]/[z + b] — ► e(— (a — b)/2u{) as Im(z/u;i) — ► oo. This 
implies in f)3. 1 j) 



ft ^r - eC-c/a*), I^i±i -> eC-c/2^) (i = 1, . . . ,m), 

jj[ [Z ~ Xj\ [Z - Xi\ 

as Im(z/ui) — > oo. Hence we obtain (|3.3|) from (|3.1|) . In the rational case [z] = z, formula 
is derived by a simple limiting procedure z — > oo. 
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3.2 Case of type A 

We apply the key identities (|3.2p and (|3.3p for studying kernel functions. For two sets of variables 
x = (xi, . . . , x m ) and y = (y±, . . . , y n ), we consider the following meromorphic function: 

jj{ [z-vj] fJi [z-yi + v] 

where k, X and v are complex parameters. Then by Proposition 13.11 this function F(z) is 
expanded as 

yiz-x^ j-, [x -x J+K ] ^ [x -y l + v + X] 

i [A]y^ ^~ j/fc + t, + ^ TT ^ - - + A ^ TT ~ x j ~ v + K 1 

fc=l L KKn;Uk Lw WJ 7=1 Ltf J J 



^ [xi-Vi + v] 

z-yk + v + c] -fj [y k - xj - v + k] 



f^i [z-Vk + v] f = A [y fc 



x,- — f 



where c = mn+nX, with the coefficients ^4i(x; k), Ak(y; X) of Ruijsenaars operators in x variables 
and y variables. By Lemma [3~T1 under the balancing condition c = run + nX = 0, we have 



Mg>(*; lXi . n+v] + wX>fo A >n R _,._„] - »• 

Also, when [u] is a constant multiple of sin(7ra/wi) or u, we always have 

H f>(x; «)n [ ","" + : + i A1 + wix* A)n ^-^-" + ;i = ^ + ^. 

i=i /=i [Xi + ^ k=i j=i [yfc UJ 

We now try to find a function y) satisfying the system of first-order difference equations 
T^Hx; y) = f[ [x ]- yi + V + X] Hx; y) (i = l,...,m), 



T^(x; y) = f[ [Vk V + K] *(z; y) (k = l,...,n), (3.4) 

fj{ Wk ~ Xj - v\ 



where r stands for the unit scale of difference operators for y variables. In order to fix the idea, 
assume that the balancing condition mn + nX = is satisfied. Then, any solution of this system 
should satisfy the functional equation 



Til 

i=l k=l 



K ]J2Mx;K)Tl<S>(x;y) + [A]^4 fe (y; X)T^(x;y) = 0, 



namely, 

[K]Di 5 ' K ^(x;y) + [X}D^ x ^( X ;y) = 0. 
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The compatibility condition for the system (|3.4p of difference equations is given by 

[xj - yk ~ t + v + A] \xj -yk + v] = [Vk ~ Xj - 5 - v + k] [y k - Xj-v) 
[xi - yk - t + v] [xi - yk + v + A] [y k - Xi - 5 - v] [yk - Xi - v + n] 

for i = 1, . . . , m and k = 1, . . . , n. From this we see there are (at least) two cases where the 
difference equation (|3.4|) becomes compatible: 

(case 1) : r = —5, A = —k, v : arbitrary, 
(case 2) : r = k, A = 5, v : arbitrary. 

In the first case, the difference equation to be solved is: 

[xi-yi+v- k] 



!/) = II r 1 i y ) (i = l,...,m), 

m r 



[Vk-Xj-V + K] 

$>{x;y) (fc = l,...,n). 



Jk - Xj - v\ 
This system is solved by 

Hence we see that $^(x;y|5, k) satisfies the functional equation 
l/IM- ^-^^(x; y\6, k) = 0, 

under the balancing condition (m — n)K = 0, namely, m = n. By setting <&a{x; y) = (x; —y), 
we obtain 

£>tt">$ A (x ;i ,|a,«) - K ^ A (x; y\6, «) = 
as in Theorem 12.11 The second case 

z£*(*; y) = n [ T" ;+ I + / ] ^ ; y ) (i = i,...,m), 

^(g;y)=n ^" g, '" t,+ i ,6l ^y) (* = !>•..,») 
y fj[ [yfc - Xj - v\ 

is solved by 

m n 

^A{x;y) = YUJ[xj-yi+v]. 

Hence, under the balancing condition run + n5 = 0, we have 

[K]D^ A {x-y) + [5}D^ A (x-y) = 0. 

This completes the proof of Theorem 12.11 

When [u] is a constant multiple of sin(7ru/o;i) or u, for any solution $(x; y) of the system (I3.4D 
of difference equations we have 



m 

i=l fc=l 



]^A l (x ] K)T^(x-y) + [\]^A k (y-,\)T« k <S>(x-,y) = [tuk + nA]$(x; y), 
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namely, 

[K]Di S ' K) <f>(x;y) + [X]D^ T ' X ^(x;y) = [ran + n\]<f>(x;y) 

without imposing the balancing condition. This implies that the functions & A (x;y\5, n) and 
^ A (x;y) in these trigonometric and rational cases satisfy the functional equations 

Di 5 > K ^ A (x; y\5, «) - D^1> A (x; y\S, k) = [(m ~ n)/C U A (x; y\S, K ) 

[K\ 

and 

{K]Di 5 > K ^ A (x;y) + [5}D^ 5 ^ A (x;y) = [m« + nS)^f A (x; y). 
respectively, as stated in Theorem 12.21 

3.3 Case of type BC 

In this BC case, we assume that [u] is a constant multiple of one of the following functions: 

u (p = 1), sin(iru/uji) (p = 2), e (au 2 ) a(u; $7) (p = 4). 
In order to discuss difference operators of type BC, we consider the meromorphic function 

2p 

s=l 

By Proposition 13.11 it can be expanded as 



m r I I rn r i 



i=l 1 



] 4 [z + xt] 



y^ [z-y k + v + c] + ^ [z + y k + v + c] 
h Iz-Vk + v] ^ k t[ [z + y k + v] ^ 

into partial fractions, where 

2p p 

c = Iijik + 2nA + E^ s ~ § + «) + E^ s = 2mK + 2n A + c (At|5 ' K) . 

S=l 8=1 

Also by Lemma |3. 11 we see that expression 

E E o* + E R r+ E * 

l<i<m;e=± l<fc<n; e=± 1<>*<P 1<^<P 

reduces to when the balancing condition c = is satisfied, or to [c] when [u] is trigonometric 
or rational. A remarkable fact is that, if the parameter v is chosen appropriately, then the 
expansion coefficients of F(z) are expressed in terms of the coefficients of Ruijsenaars operators 
of type BC. 
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Proposition 3.2. Whenv = ^(5 — X), the function F {z) defined by (|3.5p is expressed as follows 
in terms of the coefficients of Ruijsenaars operators: 

[c] n [z + Vs] m r, . , n r , 

1/, Mr 1/ \i \z±x j \ 11 [«±W + w] 

+ wE"(W) ' z i!j ( f f r r itr l A ° fe " | ^ ) 

r=1 + 2 (O — LO r )\ 

+ [ A ]2^ e (2 c ^) r w „ ^(y^k, A), 



where 



c = 2mK + 2nX + c^ 5 ' K \ t = k + \-5, v = {m — — v). 



Proof. The expansion coefficients in (|3.6h are determined from the residues at the correspond- 
ing poles. We first remark that 

[ K ] EI i x i + Ms] ro i m r , , n r , xl 

p+ _ s=i iggj + kJ rr Fj ± gj + «j tt R ± 2/* + v + A] 

* Ar 1/r Mr 1/ \1 \2Xi] *~ *~ \Xi ± X,'l ^ [sC* =t «/ + ul 

I! [Xi + ±(£-W B )][Si + i(K- W a )] J# JJ '=1 1 

8=1 

Since 

[2Xj + k] _ t4 [xj + |(w - CJ S )] 
[2xi] " A = l [ Xi - i w J ' 

we obtain 

[«] EI R + Ms] m r -I n r . i \ 1 

p + _ s=i tt R ± Xj + k\ yr [Xj ± y t + v + AJ 

i rrr 1 k , i /s \Ah [xi±xj] f\ [xi±yi + v] 

8=1 

Note that P~ is obtained from by replacing Xj with — xj (j = 1, . . . , m). We next look at 
the coefficient : 

w n bfc + Ms - v] 

Q+ = 8=1 [ggfc + AJ 

n k + i(« - w.) - fib/* + §(« - ^) - [2yfc] 

s=l 
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x TT [Vk ± Vl + A] yr [y k ±Xj+K-v] 

In view of 

[2y k + X] = ^ [y k + ±(\-u Js )} 

[2yd [y k - \u s ] 

we set — v = |A, namely, ?; = ^(S — A). Then we have 

2p 

[A] EI bfe + - v] , m r 



n + _ s=l tt [Vk ; ffl ~r A J "|~f 

^ ~ P . . 11 U„ + ,„1 11 



s=l 



v = ~{8 — A), r = k + A - <5, u = (fj,t — v, . . . , /j,2 P — v). 

The coefficient Q^j" is obtained from Q^" by the sign change yi — > —yi- The coefficient R r is 
given by 

2p i 

R Bt_ ^ TT U± ZM X A + "] fT [h(<*r-6)±Vl+V + A] 

" n[i«r.in[^. + i(*-«)]M SK-^**] fj [*(«*-*)±tt+«] 

where u; rs = w r — u; s . When u = i(<5 — A), we have 

(ui. -4- A") + «/l 

e(A?? r ) 



[±(w r - 5) ± y z + u + A] _ [±(w r + A) ± y { ] 



[i(w r -5) ±y t +v] [\{u r - A) ±y t ] 
by the quasi-periodicity of [u]. In this way y variables disappear from R r : 

n n [fa* + i(« - <s)]i=i [ ^ (u;r " 5) ± Xj] 

Note that the exponential factor e(n\r] r ) is nontrivial only in the elliptic case. In any case, from 
2n\ = c — 2rriK — c^ S ' K ' we have 

R r = [K]e(^cr] r )A^,(x; fi\5, k). 

Finally, when v = ^(5 — A) and r = k + A — 5, we can rewrite SV as 

5 ^[2 K - ") + m ^ + w) ± n ( i ( ^ _ w) ± yt + ^ + A] 

sjtr s=l 
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8 

e(mKT] r ) n [\{u) r - t) + n s - v] n n, s , xl 
_ s=l TT [gA^r — T) ±yi-t A\ 

rrri i A ri , U\ \i f-i [\( u r -t)±Ui) 

s^r s=l 

Note that x variables have disappeared again by the quasi-periodicity of [u]. Since 2mn = 
c - 2n\ - M 5 ^ and M 5 > K "> = C H T ' A ), we have 

S r = [\]e{\cq r )A° r (y;v\T, A). 

This completes the proof of proposition. ■ 

In what follows we set v = \{5 — A) and r = re + A — 5. Then Proposition 13.21 combined with 
Lemma |3.1[ implies 

m Ae, \ X \Ty l ex i ± yi + h( s + x )] i r iV^O/ IJt \ 

« z_> ^( aj ;/4 y »«)llr — z — nr^ — rz; + m> A Ax; k) 

l<r<m; e=± i=l L 1 31 2 V /J s =l 

+ [A] £ ^tte 1r, A)fi t + ■ + III + [A '^^ (K y|T - A) 

l<fc<n; e=± J=l ± X J + 2^ ~ K JJ r =l 

= 0, (resp. = [2mK + 2nA + c (/i|5 ' K) ]), 

when the balancing condition 2mK + 2nA + c^' <5,K - ) = is satisfied (resp. when [u] is trigonometric 
or rational.) Hence we have 

Proposition 3.3. Suppose that the parameters 5, k, t, A satisfy the relation S + r = re + A, and 
define v = {v\,.. .,v 2p ) by 

v s = Us - \{f> ~ A) = n s + \{t - re) (s = 1, . . . , 2p). 

Let $(x;y) any meromorphic function in the variables x = (xi, . . . ,x m ) and y = . . . ,y n ) 
satisfying the system of first- order difference equations 

T Xi $(x;y) = l[ r — -frz rr7$(z;y) (» = l,...,m), 

f^ l [x i ±y l + ^(5 - \)\ 

(1) // i/te balancing condition 2mK + 2?iA + c^'' 5 '^^ = is satisfied, then $>(x;y) satisfies the 
functional equation 

[ K ]E^ s ^(x;y) + [X]E^<S>(x;y) = 0. 

(2) If [u] is a constant multiple ofsm(iru/uJi) oru, then <&{x\y) satisfies the functional equation 

[ K ]E^ S ' K ^(x;y) + [X]E^ T ' x ^(x;y) = [2mre + 2n\ + c M5 ' k) ]<Z>(x; y). 

In fact there are essentially two cases where the system of first-order linear difference equa- 
tions (|3.7p become compatible: 



(case 1): r = —6, A = — k, u s = p s — ^(5 + re) (s = 1, . . . , 2p); 

(case 2): r = re, A = 5, u s = p s (s = 1, . . . , 2p). (3.8) 
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In the first case, the system (j3.7|) of difference equations to be solved is: 

<*-> »>• 

It is solved either by the function 
or by 

m n 

$Bc{x;y\5,K) = ~[[ll 11 G{e lXj + e 2 yi + \{5- k)\S). 

j=l/=lei,e2=± 

Hence we see that, under the balancing condition 2(m — n)K + c^' <5,K - ) = 0, $ B c(%\ y\S, k) satisfies 
the functional equation 

E^<5> BC (x; y\5, k) - E^~ K ^ BC (x; y\5, k) = 0, 

for u s = /i s — g {$ + k) (s = 1, . . . , 2p). By symmetry (|2.12p with respect to the sign change of 
parameters, we obtain 

E^^ BC {x- y\5, k) - E^^ BC {x- y\S, k) = 0, 

for v s = 2 + k) — p s (s = 1, . . . , 2p). The second case 

Zj*(z; y) = fl ^f^^ Hx; y) (i = l,...,m\ 
iJi t Xi ± yi * 

T vl^ y) = n ^flCi" 1 ^ y) (A; = 1; ■ ■ ■ ' n) 

is solved by 

m n 

^Bc(x;y) =]JY[[x j ±y l ]. 

j=ll=l 

Hence, under the balancing condition 2ttik + 2nd + c^ S ' K ^ = 0, this function satisfies the func- 
tional equation 

M^ w fflc(i;») + = 0. 

This completes the proof of Theorem 12.31 When [u] is trigonometric or rational, Proposi- 
tion [231 (2) implies that these functions $bc(x; y\5, k) and ^Bc( x ]y) satisfy the functional 
equations as stated in Theorem 12.41 
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3.4 Difference operators of Koornwinder type 

In the rest of this section, we confine ourselves to the trigonometric and rational BC cases and 
suppose that [u] is a constant multiple of sin(7ru/a;i) or u. We rewrite our results on kernel 
functions for these cases, in terms of difference operator 

m m 

D MS, K ) = J2Af(x^\S,K){Tl - 1) + J2A7( X ;v\6,k)(T- 5 - 1) 

i=l i=l 

of Koornwinder type. We remark that this operator has symmetry 
dM-S,-k) _ t~)(-H<5,k) 

X X 

with respect to the sign change, as in the case of E^ S ■ We show first that D x K differs 
from Ex K ^ only by an additive constant in the Oth order term. 

Lemma 3.2. 

(1) The constant function 1 is an eigenfunction of E± : 

(2) The two difference operators Dx K ^ and E x K ^ are related as 

E (p,\5,n) = D W, K ) + c (ji\S,k) + 1 Q 2mK + c (m|5,«)] _ [ c (mI«/0]). 

Proof. Since D x ' = £4^' — E^ S (1), statement (2) follows from statement (1). For the 
proof of (1), we make use of our Theorem 12.41 This theorem is valid even in the case where the 
dimension m or n reduces to zero. When n = 0, Theorem 12. 4| (1) implies 



[n]E^ s ' K \l) - [k)C {u ^ k) = [2mK + 



.(mI<5>«)i 



with v s = i(<5 + k) — ia s (s = 1, . . . , 2p), since &bc{x] y) hi this case is the constant function 1. 
Also from the case m = n = we have 

^ C W,k) _ [ K p(»\S,K) = [ C (A*I«,«)]. 

Combining these two formulas we obtain 

[k]E^ S ' k \i) = [k]C^ S ' k) + [2mn + c M5 ' K) ] - [ c 0*l*.*)]. ■ 

Let us rewrite the functional equations in Proposition 13.31 in terms of the operator D X ^ S ' K \ 
In the notation of Proposition 13.31 (2) we have 

[ K ]EM 5 >«)$( X ;y) + [X}E^ T ' X ^{x;y) = [2ms + 2nA + cM s >")]$( X ; y). (3.9) 

As special cases where (m,n) = (m, 0), (0,n) and (0,0), we have 

[/t]£*'(l) + [A]C Wr ' A) = [2itik + c ^ s '% 
[ K ]CM«.*) + [A]^l r ' A )(l) = [2nA + C WW], 
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and hence 

-[K]E^ K \l)^{x;y) - [A]C Mr ' A) $(x;2/) = -[2mn + c^ 5 ^]${x; y), 

[ K ]C^ 5 ' K ^(x; y) + [\}C^ T > x) $>{x- y) = [c ( ^' k) ]$(x; y), (3.10) 
By taking the sum of the four formulas in (|3.9|) and (|3.1U|) . we obtain 
[K]D^ S ' K ^(x;y) + [X}D^$(x;y) = C<S>(x;y), 

C = [2rriK + 2nX + c M<5 ' k) ] - \2niK + c if1 ^} - [2nX + c^ S '^} + [c MS '% 

In the rational case, it is clear that C = 0. In the trigonometric case, this constant C can 
be factorized. In fact, if we choose [u] = 2\J — 1 sm(iru/uji) = e(u/2uj\) — e(—u/2oji), we have 
a simple expression 

C = [mK][n\][mK + nX + M s >*\ 

Hence, under the assumption of Proposition 13.31 (2), we have 

[ K ]D^ S ' K) ^{x;y) + [X]D^ T ' X) <^(x;y) = [mK][nX][mn + nX + c W ' k) ]$(x; y) (resp. = 0), 

when [u] = e{u/2u}\) — e(—u/2ui) (resp. when [u] = u). Applying this to the two cases of (|3.8j) . 
one can easily derive functional equations as in Theorem 12.51 

4 Kernel functions for g-difference operators 

In the trigonometric case, it is also important to consider g-difference operators passing to 
multiplicative variables. Assuming that [u] = 2y/—l sin(-7rn/a;i) = e{u/2uj\) — e(—u/2u)\), we 
define by z = e(u/ui) the multiplicative variable associated with the additive variable u. When 
we write [u] = zz —z~z = — z~? (1 — z), we regard the square root z? as a multiplicative notation 
for e{u/2oj\). We set q = e(5/u>i) and t = e(K/cJi), assuming that Im(<5/o»i) > 0, namely, |g| < 1. 

4.1 Kernel functions for Macdonald operators 

We introduce the multiplicative variables z = (z±, . . . , z m ) and w = (w\, . . . , w n ) corresponding 
to the additive variables x = (x±, . . . , x m ) and y = (yi, . . . , y n ), by Z{ = e{xi/uj{) (i = 1, . . . ,m) 
and Wk = e{yk/u)\) (k = 1, . . . , n), respectively. 

In this convention of the trigonometric case, the Ruijsenaars difference operator Dx' K ^ of 
type A is a constant multiple of the (/-difference operator of Macdonald: 

m 

D (s, K)=t -i {m -i )v ( q ,t) j p[^) = ]T H t -^i TqtZi , 

i=ll<j'<m; j^i 1 ^ 

where T qZi denotes the g-shift operator with respect to Zf 

T q , z J(zi, ...,z m ) = f(zi, . . .,qzi, ...,z m ) (i = l,...,m). 

In what follows, we use the gamma function G-{u\8) of Section[2l (|2.3|) . Then our kernel function 
&a{%] y\5, k) with parameter v = k is expressed as follows in terms of multiplicative variables: 

<S> A (x; V\S, k) = e{^{f))( Zl ■ ■ ■ z m )25(w 1 ■ ■ ■ w n )^U(z; w\q, t), 
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m n 



The functional equation of Theorem 12.21 (1) thus implies 

i ira-n 

V^U(z; w\q, t) - t m - n V^U{z; w\q, t) = U(z; w\q, t), (4.2) 
which can also be proved by the expansion formula 

U{z;w\q,t)= b x (q,t)P x (z\q,t)P x {w\q,t) (4.3) 

Z(A)<min{m,ra} 

of Cauchy type for Macdonald polynomials p3]. (Formula (j4.2|) already implies that H(z; w\q,t) 
has an expansion of this form, apart from the problem of determining the coefficients b\(q,t).) 
On the other hand, kernel function ^^(x;y) with parameter v = is expressed as 

m n 

— — — — TTTT 

y A {x;y) = (zi ■ ■ ■ z m ) 2( Wl ---w n ) 2| [(zj-wi). 

j=U=l 

Then the functional equation of Theorem 12.21 (2) implies 

m n 

( ( i _ t)V toA _ (i _ q)V M _ ( i _ iY))nn^ - = o. 

3=11=1 

This formula corresponds to the dual Cauchy formula 

m n 

titled ~ w l)= E {-^P\{z\q,t)P x .(w\t iq ), 

3=11=1 Ac(ra m ) 

where A* = (m — \' n , . . . , n — A^) is the partition representing the complement of A in the mx n 
rectangle. 

These kernel functions for Macdonald operators have been applied to the studies of raising 
and lowering operators (Kirillov-Noumi [U[B], Kajihara-Noumi [S]) and integral representation 
(Mimachi — Noumi [16], for instance). We also remark that, in this A type case, a kernel function 
of Cauchy type for g-Dunkl operators has been constructed by Mimachi — Noumi |17j . 

4.2 Kernel functions for Koornwinder operators 

We now consider the trigonometric BC case. Instead of the additive parameters (fj,\, //2, A*3, Hi) 
(p = 2), we use the multiplicative parameters 

(a,b,c,d) = (e(//i/o;i),e(^2/wi),e(jU3/a;i),e(//4/a;i)). 

These four parameters are the Askey-Wilson parameters (a, b, c, d) for the Koornwinder polyno- 
mials P\{z; a, b, c, d\q, t). 

In this trigonometric BC case, the difference operator 

m m 

D m*) = -£ A +(x;li\6,k)(tI - 1) + £\4r (x; »\5, k) (T~ s - 1), 

1=1 i=l 
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which we have discussed in Section 13.41 is a constant multiple of Koornwinder's g-difference 
operator [TU]. Let us consider the Koornwinder operator 

m m 
V (a, b ,c, d \ q ,t) = J2Af(z)(T q , Zi - 1) + J^A" Wte " 1) 
i=l i=l 

in the multiplicative variables, where the coefficients Af(z) = Af(z;a,b,c,d\q,t) are given by 



a+i z \ _ (1 ~ azjjjl - bzj)(l - czj)(l - dzj) -p-r 1 - tz{Z 



3 



1 - Zizf 1 



(abcdq-l)h™-i(l - zf)(l - qzf) ]+i - 

and A~{z) = Afiz- 1 ) for i = 1, . . . ,m. Then we have D { x ^ r) = _£>( a Ac,%,^ Note that thig 
operator x>[ a ' b ' c ' d ^ q ' t ^ is renormalized by dividing the one used in [10] by the factor {abcdq~ l )^t m ~ l . 
In what follows, we simply suppress the dependence on the parameters (a,b,c,d\q,t) as V z = 
jj(a,b,c,d\q,t) ^ w j ien we re f er ^ Q operators or functions associated with these standard parameters. 

In Section [2j we described two types of kernel functions (12.141) and (I2.15|) of Cauchy type. 
Also, depending on the choice of G(u\5) we obtain several kernel functions for each type. From 
the gamma functions Gzf(u\5) of (|2.3j) . we obtain two kernel functions of type (12. : in the 
multiplicative variables, 

m n I 1,1 ±1 \ 

/p-r-r-i-r [Q 2 t?ZjW, :q) 

*„(*; w\ g , t)=( Z1 ... ^nn , i i ± r . ^ 

j= U =1 {q*t 2 Zj w t ;q) r 



and 

{q^t^z~ l wf l ;q) 
fMAiqn-h^w^^q) 



w\q, t) = ( Zl ■ ■ ■ *mr n/? ni L l,-liZ\.±^ 



respectively, where we put (3 = k/5 so that t = q@. Similarly, we obtain two kernel functions of 
type (l2~TSl) from G±(u\5): 

m n 

* + (z- t w\ q ,t) = e (^)nn n 

i = li=lei,e2=± 
m n 

M^k,i)=e(-snn n ^m^t^ 

j = li=lei,£2=± 

m n 

where f(x;y) = n ^ x 2 + m ^ yf + ^(k 2 — <5 2 ). Each of these four functions differs from 
j=i i=i 

the others by multiplicative factors which are 5-periodic in all the x variables and y variables, 
It should be noted, however, that they have different analytic properties. In the following we 
denote simply by &(z;w\q,t) one of these functions. 

The kernel function ^(z; w) = $?bc( x ] y) of dual Cauchy type is given by 

m n m n 

*(z; w) = + zj 1 - Wl - wf 1 ) = Ullizj - Wl ) (1 - zfw' 1 ) , ( 4 -5) 

3=11=1 j=ll=l 



which is precisely the kernel function introduced by Mimachi |15j . 

For the passage from additive variables to multiplicative variables, we introduce the multi- 
plicative notation for the function [u]: For z = e(u/u}\), we write {z) = [u]. Namely, we set 



i i 



(z) = z 2 — z 2 = —z 2 (1 — z ), z = e(u/uji), 
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with the square root regarded as the multiplicative notation for e(u/2ui). This function (z) 
is a natural object to be used in the case of BC type, because of the symmetry (z^ 1 ) = —(z). 
In this notation, the coefficients Af{z) of the Koornwinder operator T> z are expressed simply as 

A + ( x _ {aZj){bZi){cZi){dZi) -i-r {tZi/Zj){tZiZj) _ 

' )= WW) , JL* <«/*>)<**>) ( - 1 -" m >- 

4 

It should be noted also that our parameter c^' 5 ' K ) = ^2 fj, s — (5 + k)+loi passes to multiplicative 

8=1 

variables as 



4 

u + S 

s=l 



-(zabcd/qt), z = e(u/ui) 



with a minus sign. Then, Theorem 12.51 can be restated as follows. 
Theorem 4.1. 

(1) The function &(z;w\q,t) defined as above satisfies the functional equation 

(t)V z <S>(z;w\q,t) - (t)V w <S>(z;w\q,t) = it m )(r n )(abcdq~ l t m ~ n ~ l )^{z]w\q,t), 

where T> w denotes the Koornwinder operator in w variables with parameters (a, b, c, d) 
replaced by (y/qt/a, y/qt/b, y/qt/c, %fqb/d). 

(2) The function \V(z;w) defined as (|4.5p satisfies the functional equation 

{t)V z *(z;w) + (q)V w V(z;w) = {t m ) (q n ) {abcdt™- 1 q n - x )V> {z; w), 

where D w denotes the Koornwinder operator in w variables with parameters (a, b, c, d\t, q). 

Statement (2) of Theorem 14.11 recovers the key lemma of Mimachi [151 Lemma 3.2], from 
which he established the dual Cauchy formula 

m n 

YlYl(zj + zr 1 - wi - W7 1 ) = Yl {-l) lX * l P\(z;a,b,c,d\q,t)P x *(w;a,b,c,d\t,q) 
3=11=1 Ac(n m ) 

for Koornwinder polynomials, where the summation is taken over all partitions A = (Ai, . . . , A m ) 
contained in the m x n rectangle, and A* = (m — X' n , . . . ,m — A^). By this formula, he also 
constructed an integral representation of Selberg type for Koornwinder polynomials attached 
to rectangles (n m ) (n = 0,1,2,...). We expect that our kernel function $(z;w\q,t) of Cauchy 
type could be applied as well to the study of eigenfunctions of the (/-difference operators of 
Koornwinder. As a first step of such applications, in Section [5] we construct explicit formulas 
for Koornwinder polynomials attached to single columns and single rows. 



5 Application to Koornwinder polynomials 

In this section, we apply our results on the kernel functions for Koornwinder operators to 
the study of Koornwinder polynomials. In particular, we present new explicit formulas for 
Koornwinder polynomials attached to single columns and single rows. 

To be more precise, we make use of the kernel functions to express Koornwinder polyno- 
mials Pnr\(z; a, b, c, d\q, t) (r = 0, 1, . . . , m) and Pm(z; a, b, c, d\q, t) (I = 0, 1, 2, . . .) in terms of 
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certain explicitly defined Laurent polynomials E r (z;a\t) and Hi(z; a\q,t), respectively (Theo- 
rems ED and • We remark that these Laurent polynomials E r (z;a\t) are Hi(z;a\q,t) are in 
fact constant multiples of the BC m interpolation polynomials of Okounkov [18j attached to the 
partitions (l r ) and (I), respectively. (This fact will be proved in AppendixO) Namely, Theo- 
rems 15.11 and 15.21 provide with two special cases of the binomial expansion of the Koornwinder 
polynomials in terms of BC m interpolation polynomials as is discussed in Okounkov |18) and 
Rains [19] . 

Once we establish the fact that E r (z;a\t) and Hi(z;a\q,t) are interpolation polynomials, 
Theorems 15.11 and 15.21 can also be obtained from Okounkov's binomial formula [18] , together 
with Rains' explicit evaluation of the binomial coefficients for the cases of (l r ) and (I) [190. 

Before starting the discussion of Koornwinder polynomials, we introduce a notation 

(z; w) = (zw)(z/w) = z + z~ x — w — w -1 

which corresponds to [it±«] = [u+v][it — v] in additive variables. This expression, which appears 
frequently in the discussion of type BC, deserves a special attention. Note that 

(a -1 ; 6) = (a;6 _1 ) = (a;b), (b; a) = - (a; b) , (a; b) + (b; c) = (a; c) , 

as clearly seen by the definition. Also, the Riemann relation for [u] can be written as 

(z;a) {w;a) (z;w)(a;b) 



(z; a) (b; c) + (z; b) (c; a) + (z; c) (a; b) = 0, 



(z;b) (w;b) (z;b)(w;b}' 



5.1 Koornwinder polynomials 

We briefly recall some basic facts about Koornwinder polynomials; for details, see Stokman [23] 
for example. 

Let K = Q(a2, 62,02,(^2,52^2 J be the field of rational functions in indeterminates, rep- 
resenting the square roots of the parameters a, b, c, d, q, t, with coefficients in Q, and 
Kfz^ 1 ] = IfCfz^ 1 , . . . , z^\ the ring of Laurent polynomials in m variables z = (z±, . . . , z m ) with 
coefficients in K. Then the Weyl group W = {±l} m x 6 m of type BC (hyperoctahedral group) 
acts naturally on K.[z ±1 ] through the permutation of indices for the z variables and the individual 
inversion of variables Z{ (i = 1, . . . , m). The Koornwinder polynomial P\(z) = P\(z; a, b, c, d\q, t) 
attached to a partition A = (Ai, . . . , A m ) is then characterized as a unique VF-invariant Laurent 
polynomial in K^ 1111 ] satisfying the following two conditions: 

(1) P\{z) is expressed in terms of orbit sums m^z) = Yl z u as 

veW.fi 

P\(z) = m\(z) + y^cx^m^z) (c A)A , G IK), 

where < stands for the dominance ordering of partitions. 

(2) P\(z) is an eigenfunction of Koornwinder 's (/-difference operator D z : 

V z P x (z) = d x P x (z) for some d x £ X. 



x The authors thank Professor Eric Rains for pointing out this connection with the interpolation polynomials 
and the binomial formula. 
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These polynomials P\(z), indexed by partitions A, form a IK-basis of the ring of VF-invariants 
[z ±1 ] M/ . Also, the eigenvalues d\ are given by 



-1 j.— m+i\ 

i=l i=l i=l 



in our notation {z;w) = z + z 1 — w — w 1 , where a = (abcdq 1 ) 2 . (Note that Po(z) = 1, and 
d = 0.) 

We give below new explicit formulas for Koornwinder polynomials P^r^z) attached to sin- 
gle columns (l r ) (r = 0,1,..., m), and P^(z) attached to single rows (r = 0,1,2...). As 
we already mentioned, our explicit formulas provide with the two special cases of Okounkov's 
binomial expansion of the Koornwinder polynomials in terms of BC m interpolation polyno- 
mials. We also remark that, in the cases of type B, C, D, some conjectures have been proposed 
by Lassalle |11| on explicit formulas for Macdonald polynomials attached to single rows. The 
relationship between his conjectures and our approach will be discussed in a separate paper. 

In order to formulate our results, we define a set of VT-invariant Laurent polynomials E r (z; a) 
with reference point a (r = 0, 1, . . . , m) by 

E r (z;a\t)= ]T (z ll ;fi- 1 a)(z i2 ;f*~ 2 a)---(z lr ;fr- r a). (5.1) 

l<il <---<i r <m 

As we will see below, these Laurent polynomials are VF-invariant in spite of their appearance, 
and they can be considered as a variation of the orbit sums m^r^z) (r = 0, 1, . . . , m) attached 
to the fundamental weights. In fact, these Laurent polynomials E r (z; a\t) (r = 0, 1, . . . , m) are 
essentially the BC m interpolation polynomials of Okounkov attached to single columns (l r ) (for 
a proof, see Appendix [C|). We remark that these polynomials had appeared already in the work 
of van Diejen [3] in relation to the eigenvalues of his commuting family g-difference operators 
for this BC m case. They are also used effectively by a recent work of Aomoto-Ito [T] in their 
study of Jackson integrals of type BC. 

Theorem 5.1. The Koornwinder polynomials Pn r \[z;a,b,c,d\q,t) attached to columns (l r ) 
(r = 0, 1, . . . , m) are expressed as follows in terms of E s (z; a\t) (s = 0, 1, . . . , m): 

P ( h A\ A V" (t m - r+1 it m - r *b, t m ~ r ac, t m ~ r ad) til 

P(ir)(z;a,b,c,d\q,t) = ^ {t^ m -r)abcd) t ■ E r ^{z; a\t) , (5.2) 



=0 



where (a) t ,i = (a)(ta) ■ ■ ■ (t l 1 a), and (01, . . . ,a r ) t ,i = (ai)t,l ■ ■ ■ (a r )t,i- 

1 ( l \ 1 

By using (a)t,i = (— i) t~z \2/a _ 2 (a;t)i, formula (|5.2[) can be rewritten as follows in terms of 
ordinary i-shifted factorials of [5]: 

s-^\t T ,t ab,t ac,t ad,t)i . , 
P(ir){z;a,b,c,d\q,t) = ^ -j- E r ^{z;a\t). 

1=0 tUJ + ( m ~ r V(t, t 2 ( m - r )abcd; t)i 

For the description of Koornwinder polynomials attached to single rows, we introduce a se- 
quence of W-invariant Laurent polynomials Hi(z; a\q,t) (I = 0, 1, 2, . . .) as follows: 

Hi(z;a\q,t) 

E f\ q,Vl '''fr m ^ a ^ ^ t( ? la )w ■ ■ ■ ^ t m - 1 ^ + - +u ^a) q , Um , (5.3) 

u 1 +.-+u m =l {q)q ' Vl '" [q)q ' ,/rn 
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where 

(z',a)q,l = (z;a)(z;qa) • • • (z;q l ~ l a) = (-l) l q~^) a ~ l (az,az~ 1 ;q) r 

Note that (z;a) q i is a monic Laurent polynomial in z of degree /. (The VF-invariance of 
H[(z;a\q,t) will be proved in Lemma 15.41 below.) These Laurent polynomials Hi{z;a\q,t) can 
be regarded as a BC m analogue of the A m _i Macdonald polynomials attached to single rows: 

"01 ' ,2™ . 



, 1+ ... + , m=r (^k---(^)^ 1 

Also, they are special cases of -BC m interpolation polynomials attached to single rows (I) (see 
Appendix |C|) . 

Theorem 5.2. The Koornwinder polynomials Pr r \(z;a,b,c,d\q,t) attached to rows (r) (r = 
0, 1, 2, . . .) are expressed as follows in terms of Hi(z; a\q, t) (I = 0, 1, 2, . . .): 



(t) q , r . ,, , (t^^ab^ac^ ad)^ 

^-P^aAc^t) = {q ^ m . 1)ahcdqr 



r 



X 



\- (-lYiq-^t^-Vabcdgr- 1 )^ 



Let us denote by p r (z; a, b, c, d\q) the Koornwinder polynomial Pr r )(z; a, b, c, d\q, t) in the one 
variable case (r = 0, 1,2, . . .). Note that, when m = 1, if;(z; a|q, t) reduces to {t) q! i(z;a) q j/(q) qj i. 
Hence Theorem 15.21 for m = 1 implies 

{ab,ac,ad) qtr J r (-l) l (q- r ,abcdq r - 1 ) q j 

p r {z; a, b, c, <% = — > , (z; a)^ 

(abcdif /g.rf"^ (q,ab,ac,ad)qj 



{ab,ac,ad;q) r 
a r (abcdq r ~ 1 ; q) r 



,abcdq r ,az,a/z 
ab, ac, ad 



which recovers the well-known 4^3 representation of the (monic) Askey- Wilson polynomials. 

We prove these Theorems 15.11 and 15.21 in Subsections 15.21 and 15.31 by means of the kernel 
functions of dual Cauchy type, and of Cauchy type, respectively. 



5.2 Case of a single column 

We first explain some properties of the elementary Laurent polynomials E r (z;a\t) (r = 0,1, 
. .. ,m). 

Lemma 5.1. The Laurent polynomials E r (z;a\t) (r = 0,1,..., m) are characterized as the 
expansion coefficients in 

m m 

~[(w;zj) = y^(-l) r E T (z; a\t)(w; a) t ,m-r, 

j=l r=0 

where (w;a)t,i = (w;a)(w;ta) ■ ■ ■ (w;t l ~ 1 a). In particular, E r (z;a\t) is W -invariant for each 
r = 0, 1, . . . , m. 

Proof. Since the uniqueness of expansion in terms of {w;a)t, r (r = 0,1, . . . ,m) is obvious, we 
show the validity of the expansion formula above. Note that E r (z; a\t) can be expressed as 

E r {z;a\t) = ^JJ(^;tl J <*la), I c <i = {j £ {l,...,m}\J \j<i}. 

\I\=rieI 
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Hence we have the recurrence formula 

E r (z;a\t) = (zi;a)E r ^i(z';a\t) + E r (z';ta\t), z' = (z 2 , ■ ■ ■ , z m ). (5.5) 

By using this recurrence, one can inductively prove the expansion formula of lemma. In fact, 
by (w; zi) = (w;a) — (zi;a), we compute 

m m m 

X\{w,Zj) = (w;a)Y[(w,Zi) - {z\\ a)\\{w; z^) 

j=l i=2 i=2 

m—1 

= (w; a) ^2 (- l Y' E r(z'; ta\t)(w; ta) t , m -i- r 

r=0 
m—1 



- (zi; a) (-l) s E s (z; a\t)(w; a) t , m -i- s 

s=0 

in 

^2(-l) r (E r (z?;ta\t) + (z 1 ;a)E r _ 1 (z';a\t))(w;> 

r=0 
m 

^2(-l) r E r (z;a\t){w; 



r=0 

We now proceed to the proof of Theorem 15.11 Mimachi's dual Cauchy formula for Koorn- 
winder polynomials can be written as 

m n 

Y{Y[(wi; zj) = ^2 {-l) lXl P\(z;a,b,c,d\q,t)P x *(w;a,b,c,d\t,q). 

3=11=1 Ac(n m ) 

When n = 1, this formula implies 

m m 

X\{w,Zj) = ^2(-l) r P {in (z\q,t)p m _ r (w\t), (5.6) 

3=1 r=0 

where we have omitted the parameters (a,b,c,d). Namely, the Koornwinder polynomials at- 
tached to single columns are determined as expansion coefficients of the kernel function for 
n = 1 in terms of the monic Askey-Wilson polynomials pi(w\t) = pi(w;a,b,c,d\t) with base t. 
On the other hand, we already have the expansion formula 

m m 

Y[(w; Zj ) = a\t)(w; a) t , m -l- (5.7) 

i=i i=o 

Recalling that 

(ab,ac,ad) t iJ^ {t~ l ,abcdt l ~ l ) t , r 
{abcdt 1 }t,if^ {t,ab,ac,ad)t, r 

E(t r+1 ,t r ab,t r ac,t r ad)ti-. ri . . , 

we consider to express {w;a) t j in terms of Askey-Wilson polynomials p r (w\t). 
Lemma 5.2. For each I = 0, 1, . . ., one has 

w *(t r+1 ,t r ab,t r ac,t r ad) f ;_ r , , . , . 

(„;„),,, = g(-l)'- ^ \ t ^ cit ,; ) J'' V MO = 0,1,2,...). (5.9) 
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We omit the proof of this lemma, since it can be derived as a special case of the connection 
formula for Askey- Wilson polynomials with different parameters (see |U). Note that, if we set 
d = t 1 ' 1 1 a in (|5.8p . then pi(w; a, b, c, t 1_l /a\t) = (w; a)t,i- 

By substituting (|5.9p into (|5.7p . we obtain 

TT/ \ ^ r ly gl + \t m - r ab,t m ~ r act™-' a^r-l „ , ... , ... 
j=i o<i<r<m \z,aocazy >) t:r _i 

Comparing this formula with (|5.6p we obtain 

r Um—r+l +m—r r ,u +m—r r ,„ +ra—r r , r i\ 

p a-)(*l**) = g )), r _ t ^ (z;a|t) ' 

as desired. 



5.3 Case of a single row 

Recall that the kernel function of Cauchy type 

a-t-T-t-r (q 2 t2ZjW, ;q) 
Mz; W \ q ,t)=( Zl ...z m )^l J * / 1 flf , 

j=ii=i(? 3 * ;?)«> 

defined in (|4.4|) . satisfies the difference equation 

(t)V z ^(z;w\q,t) - (t)V w <f>(z;w\q,t) = {^{^{abcdq-H^-^Q&wfat), 

where T> w denotes the Koornwinder operator in w variables with parameters (a, b, c, d) replaced 
by (y/qt/a, y/qt/b,y/ql/c,y/qt/d\q,t). We set hereafter 

a = \Zqi/a, b = \/qi/b, c = yfqi/c, d=y/qt/d. 

Also, for any Laurent polynomial f(z) € K.fz 1 ' 11 ], we denote by f(z) € K^ 1 * 11 ] the Laurent 
polynomial obtained from f(z) by replacing the parameters (a,b,c,d) with (a, b,c, d). 

Let us consider the special case where t = q~ k (k = 0,1,2,...). Then the kernel function 
$>(z; w\q, q~ k ) reduces to a Laurent polynomial in (z,w): 

m n m n 

<S>{z;w\q, q ~ k ) = {z 1 ---z m )- kn J{\{{q^ 1 - k h jW f l - q ) k = {~l) kmn \{\{{w V , q^' k ) Zj ) q , k . 

3=11=1 3=U=1 

Ignoring the sign factor, we set 

m n 

^ k (z;w) = Hl[(w l ;q^ 1 - k h j ) qjk . 

j=ll=l 

Note that 

(w; q^ l - k) z) q , k = (w; q^~ k+1) z) {w; q^~ k+2 h) ■ ■ ■ (w; q^ k ~^ z) 

is invariant under the inversion z — > In what follows, we analyze the case where n = 1 and 
t = q~ k (k = 0, 1, 2, . . .). In this case, the kernel function 

m 

$- k (z;w) = JJ(u;;g5 (1 " fc) 3j>* 
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is a symmetric Laurent polynomial in w of degree km. Also, this kernel function satisfies the 
functional equation 

V z $_ k (z;w)-V w $_ k (z;w) = -(t m )(abcdq~H m - 2 )<5>^ k (z-w). 

With a = (abcdq" 1 )? , this formula can be written as 

V z <5>_ k {z; w) = (V w - (at" 1 - 1 ; at^))^^ k (z; w). (5.10) 

Noting that &_ k (z;w) is a symmetric Laurent polynomial, we expand this kernel in terms of 
the monic Askey-Wilson polynomials Pi(w\q) = pi(w;a,b,c,d\q) in w with the twisted parame- 
ters, so that 

m km 

<5>- k (z;w) = J\(w;q^ l ~ k) Zj ) k = ^Gi(z)p km -i(w\q). 
3=1 1=0 

The Laurent polynomials G[(z) (0 < / < km) are uniquely determined by this expansion, and 
hence VF-invariant. 

Lemma 5.3. When t = q~ k (k = 0,1,2,...), the W -invariant Laurent polynomials G\(z) 
(0 < I < km) defined as above are eigenf unctions ofT> z : 

V z Gi(z) = (at m - l q l ; at^^G^z) (0 < I < km). 

Proof. Note that 

T> w p km _i(w\q) = (aq hm ~ l ;a)p km -i(w\q) = (at' m q~ l ;a)p km -i(w\q), 

where a = (abcdq -1 )? . Hence, by a = tqz (abcd)~ z = t/a, we obtain 

V w p km -i{w\q) = (t l ~ m q~ l /a;t/a)p km -i(w\q) = (at™' 1 q l : ; at~ l )p km _i{w\q) . 

In view of (|5.10p . we compute 

(V w - (at^-^r^pkm^iwlq) = ((at^J^at- 1 ) - (at™' 1 ; odT^p^^q) 

= (at m - 1 q l ;at m - 1 )p km _ l (w\q). 

Hence, (|5.10p implies 

km km 

YPzGl(z)Pkm-l{w\q) = J2 G l( z ){Vw - (at^-^at-^p^iwlq) 
1=0 1=0 

km 

= ^(at m - 1 q l -at m ~ 1 )G l (z)p km _ l (w\q). 
1=0 

Namely, 

V s Gi{z) = (at m - 1 q l ;at m - 1 )G l (z) (I = 0, 1, . . . , km). U 

As we have seen above, each Gi(z) (7 = 0, 1, ... , km) is an eigenfunction of T> z with precisely 
the same eigenvalue as the one for the Koornwinder polynomial P^(z) attached to the single 
row of length I. At this moment, however, we cannot conclude this Gi(z) is indeed a constant 
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multiple of the Koornwinder polynomial Pm (z) specialized to the case t = q~ k . This is because 
different partitions A may give the same eigenvalue 

m 
i=l 

under this specialization. This point will be discussed later after we determine an explicit 
formula for Gi(z). 

Since we already know the relationship between the Askey- Wilson polynomials pi(w\q) and 
the Laurent polynomials {w; a) q i, we consider to expand the kernel function $>_k(z; w) in terms 
of {w,a) q j. 

Lemma 5.4. 

(1) When t = q~ k (k = 0, 1, 2, . . .), the kernel function <£_fc(z; w) has the following expansion 
formula: 



km 



q,km—l j 



(5.11) 



i=i 



1=0 



where Hi(z) stands for H\(z\ a\q, t) defined as (|5.3j) with t = q k for I = 0,1, ... , km. 
(2) The Laurent polynomials Hi(z; a\q,t) (I = 0, 1, 2, . . .) are W -invariant. 

Proof. Statement (2) follows from the expansion formula (I5.11|) of statement (1). Since 
<&_k{z;w) is VF-invariant in the z variables, formula (|5,lip implies that H[{z;a\q,q~ k ) is W- 
invariant for k > l/m. Hence we see that Hi(z;a\q,t) itself is VF-invariant as a Laurent polyno- 
mial in K[z ± ] for each Z = 0,1,2,.... 

In the following proof of statement (1), we omit the base q, and write (w; a); = (w; a) q j. We 
first present a connection formula for the Laurent polynomials (w; a)i and (w; b)i with different 
reference point a, b: 



r=0 



(q r ab,b/a) r {w;a)i_ r , 



(-1) 



r (Or 

(q)r 



which is equivalent to the g-Saalschiitz sum [4] 
(bw,b/w;q)i 



392 



q l ,aw,a/w 
ab,q l ~ l a/b 



(ba,b/a;q)i 

Let us rewrite the formula above in the form 
l 



r=0 



Hence, 



(w;q^)z)„ = ^2(-iy 

u=0 



/a) u (w;a) 



[l—V ■ 



This implies 

(w;q^ 1 -^z 1 ) lll (w;q^ 1 -^z 2 ) 



1-12 
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/'i 



i/i=0 



Mi 



( Zl -q^ X -^/a) Vl (w; a)^-^ (w; q^ 1 ^ z 2 ) ^. 



Then in each term, we expand (w;q^ 1 ^ 02)^2 m terms of {w^q^ 1 ul a)^ 2 ^ U2 (0 < v 2 < ^2) to 
get 



/mi 1 

Ml A*2 



/'2 



Ml 




M2 






.^2_ 



(* 1;9 *(l-"l)/a) 



''1 



Y.(- 1 ) V1+V2 

v\=0 1/2=0 

x {z 2 -q^ l ~^-^-^/a) U2 {w-a)^ 2 . Ul . U2 



By repeating this procedure, we finally obtain 



~[(w;q? 



i=i 



D-D M n 

1/ j=i 



Mj 



for any [i = (m, . . . , // m ), where the sum is taken over all multi-indices v — iyx^ • • • ; 

) such 

that V{ < Hi (i = l,...,m). As a special case of this formula where fj,± = ■ ■ ■ = fi m = k, 
\n\ = km, we get the expansion formula 

m km 

<f>- k (z;w) = fj(^;g^ (1 " fc) ^')fc = J^flj (*)(«;; a) fcm _j 
i=i 1=0 

for ^_ k {z;w). Here the coefficients are determined as 



W\=l 3=1 



k 



|(i-fc)-fc(j-i)+E n 
{zf,q i<j /a) Vj 



vi+-+v m =l 3=1 W/!/j 
with t = q~ k . Replacing the parameters a by a = y/qt/a, we obtain 

m km 

<5>- k (z;w) = JI(w;^ (1-fe) 2j)ft = ^2Hi(z)(w;x/qt/a) km -i, 



3=1 



where 



m (t) E 

= E n^te;* -1 ^ (* = T fc )- 

vl +...+„ m =y=i 

We now have two expansions of the kernel function &_k(z',w): 
= y^G r (z)pkm-r{w\q) = y~]Hi(z){w; \/qt/a) q ,km- 



(5.12) 



r=0 



=0 



Also, from Lemma 15.21 we see 



1 v ^i- r {<i r+1 ,q rab ,q r a c ,q ra d)q,i-r 

(w; a) q j = 2^(-!) 71 zIm " — PrMg), 



r=0 



(<7, abcdq 2r ] 



q,l—r 
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and hence, 

/ r-r, v v^/ isi- r (q r+1 , t( f +1 /ab,tq r+1 /ac,tq r+1 /ad) q> i- r „ , 
(«,; V5t/a> ffl , = g(-l) (q^q^/abcd)^ PrHq) 

for / = 0, 1, 2, . . .. Substituting this into (|5.12|) . we obtain an expression of G r (z) in term of Hi(z) 
as follows: 

- i n km-r+l j. km—r+1 / n u j-Jsm— r+1 l nn j. n km—r+l l nr i\ 

n ( \ _ sr^ ( nr -A? /ao,tq /ac,tq /aa) q , r -i . , 

Gr(Z) "^ (_) (q^-^/abcd)^ HI{Z) 

_ A , Y _t {q l ~ r t~ m , q^t^/ab, q^t^/ac, q 1 ^ t l ~ m / ad) q ^ 

^ ' faqW-'kw-^/dbcd)^ l[z) 

(t m ,t m - l ab,t m - l ac,t m - 1 ad) q , rS ^ {-I) 1 {q~ r ,t 2{ - m ' lS > abcdq r ~ l ) q ^i 
~ (q,t^ m - 1 )abcdq r - 1 ) qi r ^ {t m , t m ~ l ab, t^ac, t^ad)^ 

From the expression obtained above, it is clear that 

H[{z) = mm (z) + terms lower than (I) with respect to <, 

\Q/q,l 

and 



G r (z) = j | g ' r m( r ) (z) + terms lower than (r) with respect to <. 

Note that (i) g , r = (Q~~ k )q,r 7^ for < r < /c. Also, we already know that each G r (z) (r = 
0, 1, . . . , km) satisfies the difference equation 

V z G r {z) = {at m - l q r -at m - l )G r {z). 

Suppose in general that a partition A = (Ai,...,A m ) satisfies the condition Aj — Aj+i < k 
(i = 1, . . . , m — 1). Since 

Ai + k < A 2 + 2k < ■ ■ ■ < X m + km 

in this case, it turns out that the eigenvalue for any partition fj, < A is distinct from d\ 
when the square roots of a, b, c, d, q are regarded as indeterminates. For such a partition A, 
the Koornwinder polynomial P\(z) = P\(z;a,b,c,d\q,t) can be specialized to t = q~ k , and 
any eigenfunction having the nontrivial leading term m\(z) must be a constant multiple of 
P\(z; a, b, c, d\q, q~ k ). This implies that, for each r with < r < k, G r (z) is a constant multiple 
of P( r )(z) specialized to t = q~ k : 

G r (z) = ^L P{r) ( z )\ (0<r<k). 

\Q/q,r 

For each r = 0, 1, 2, . . ., consider the Laurent polynomial in IKfz 1 * 11 ] defined by right-hand side 
of the explicit formula (|5.4f) of Theorem 15.21 Then the both sides of (|5.4f) are regular at t = q~ k 
(k > r), and they coincide with each other for t = q~ k (k = r, r + 1, . . .). Hence the both sides 
must be identical as rational functions in frz. This completes the proof of Theorem 15.21 
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A Remarks on higher order difference operators 

In the case of type A, an explicit commuting family of higher order difference operators, denoted 
below by Dr,x { r = l,---, m ), including Dx as a first member, has been constructed by 
Ruijsenaars [20]. He also proved in |23j that, when m = n, the kernel function of Cauchy 
type, corresponding to our §a(%] y\5, n), intertwines the whole commuting families of difference 
operators in x variable and y variables. (For the comparison of our &a{x] y\8, a) with the type A 
kernel functions of Ruijsenaars [23], see Appendix IB. 11 ) 

Fix any nonzero entire function [x] satisfying the Riemann relation as in Section [2j We 
consider a sequence of difference operators Df^' (r = 1, . . . , m) defined by 

*W - E II 'T-CT' lK- (a.i) 

7C{l,...,m}; \I\=r 1 % 3i iel 

Then, from the result of [20] and its degenerate cases, it follows that these operators D^' K 
(r = 1, . . . , m) commute with each other. In this setting the same kernel function 

* r \x \ -Pf G(xj +y t + v - k\5) 
<S> A (x;y\5,K) = I I — f — - — . . 

±L G(xj + yi + v\5) 



'r,x 



as in (|2.5p for the case m = n satisfies the difference equation 

Di 5 f<S> A (x- y\S, k) = Dffl* A (x; y\S, k) (A.2) 
for all r = 1, . . . , m. This functional equation is in fact equivalent to 

E-i-r [Xj - Xj + K] -i-r [Xj+yi + V - k] 
\xi-Xj] *-L \xi + y i + v] 

\I\=rieI;j$I 1 31 iel;l<l<m 1 1 yL 1 

sr^ tt hJk -yi + n] rr [m + Xj + v - k] 
~ 2^ 11 \y k -yi\ J-l \ y +x +v \ ' 

\K\=r keK-l^K VyK y ' J keK;l<j<m yyK J 1 

which is precisely the key identity of Kajihara-Noumi [6] Theorem 1.3], that was derived from 
the determinantal formula of Frobenius. 

In the trigonometric case, it is convenient to consider the generating function 

m 

v^\u) = = £ II v^U^ 

r=0 \I\= r it I. , ; I % 3 iel 

of the Macdonald (/-difference operators, passing to the multiplicative variables. Then from the 
eigenfunction expansion (|4.3p . it follows that the kernel function II(z; w\q, t) of (|4.ip for the 
variables z = (z±, . . . , z m ) and w = (w\, . . . , w n ) satisfies the functional equation 

(iTu^ooV^^uMzMQ^) = (t n ^;t)oo^' i) (t m n)n(z; w\q, t). 

It would be an important problem to find an elliptic extension of this formula for the case 
m ^ n. 

As for type BC, a commuting family of higher order difference operators for E^ S has 
been constructed explicitly by van Diejen [SJ [3] in the trigonometric case, and inductively by 
Komori-Hikami [9] in the elliptic case. We expect that our kernel function <&bc( x ', y\5, k) should 
intertwine the whole commuting families of higher order difference operators in x variables 
and y variables (at least under the balancing condition in the elliptic case), similarly to the A 
type case. 
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B Comparison with [23 



As we already mentioned, in his series of works [23] Ruijsenaars has constructed kernel function 
for elliptic difference operators of type (A m -i,A m -i) and (BC m , BC m ) (see also [TOj [TTJ ) . In 
this section, we clarify how our difference operators and kernel functions in the elliptic case are 
related to those of Ruijsenaars. 

In what follows, we confine ourselves to the elliptic case, and set 

p = e(u 2 /oJi), q = e(5/wi); [x] = -z~* 9(z;p), z = e(x/uj 1 ) (B.l) 

as in Section [2TT1 (2) Elliptic case, assuming that lm(u 2 /u)i) > 0, Im(<5/a>i) > 0. We remark 
that the quasi-periodicity of the function [x] is described as 

[x + u r ] = e r e(7 ]r (x + ^))[x] (r = 0,1,2,3), 

where ujq = 0, u>s = —u\ — u>2, and 

e = 1, ei = e 2 = e 3 = -l; % = Vi = 0, 7 l2 = -^, V3 = ^- 

3 1 

(We use below the index instead of 4 for oj r , e r and r/ r .) Since f\ [k^s] = — 2e(— ttj-) = —2p~z 



in this case, we have the duplication formula 

3 



2 . 
r=0 



In defining kernel functions, we use the gamma function 

G(x\S) = G + (x\S) = e(^(f ))T(pz;p,q) 

of (|2,4|) . associated with [x]. 

In the works of Ruijsenaars [23] , the two periods uj\, uj2 and the scaling constant 5 are 
parametrized as 

7T 

, , (l H , , 11- 5) X « ■ )) 

u)\ — — , u>2 — ICL+ , o — za_ 
r 

in view of the symmetry between U2 and 5 (or p and q). In terms of the i?-function 

TZ(x) = "i?(r,a+;x)" = 8(ph;p) 
defined in [231 I) (1-21)], our [x] is expressed as 



x 



2^' 



Note that lZ(—x) = 1Z(x). The elliptic gamma function 

g{x) = li G{r,a + ,a^xf = T(p^ q^ z;p, q) 
of Ruijsenaars [23\ I, (1.19)] is related to our G(x\5) by the formula 

G(x\8)=e(£- i (f))g(x + ±(u2-5)). 

The function G(x) is symmetric with respect to 0J2 and 5 (p and q), and satisfies the functional 
equation 

g(x + \6) = n(x)g(x - i<j), g(-x) = g(x)-\ 
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B.l Case of type A 

As in (|A.lj) . we consider the commuting family of difference operators -Drfi^ (r = 1, . . . ,m) of 
type A in the variables x = (x±, . . . , x m ). These operators are in fact identical to the difference 
operator A r:+ (— x) of [23l I, (2.1)], up to multiplication by constants: 



"T^i^X-i Xj ~\~ tv 



under the identification of the parameter "//' = k; we denote below this operator by A r - X = 
u A r ^ + (— x)" . In the multiplicative variables Z{ = e(xi/u>i) (i = 1, . . . , m), t = e(n/uji)^ the same 
operator is expressed as 



A - Ar(m-r) d (S,k) _ V- TT 6{tZj/ Zj]p) yr 

\I\= r iei; Hi K 11 vl ' iei 
The kernel function <& A {x\ y\S, k) of (|2.5p for the m = n case can be expressed in terms of Q(x) as 

(I m m \ \ 

S A {-x;y) = J I — 1 - 

^Gi-Xj -yi-v + k + 5 (d - w 2 )j 

nT(pczjWi/t;p,q) 
— (c = e(u/wi)). 

■ j =1 r{pcZjWi;p,q) 

The second factor 5a(— x;y) of ^a^; k) coincides with u S^(—x,y)" of [231 I) (2-6)], under 
the identification of parameters "/i" = k and = — 1> + ^(<5 + « — o^)- Since the first factor of 
$a(^; y\5, k) is an eigenfunction of T*. and T* fc with an equal eigenvalue e(— mn/2uji), from (|A.2|> 
we obtain ^l r . ! _ a; 5A(— x; y) = ^l rj _ y 5 J 4(— x; y), namely 

Ar >x S A (x; y) = Ar- y S A (x; y) (r = 1, . . . , m). 

This gives formula (2.5) in [23\ I] with ll 5 = +". The statement for "5 = — " follows from the 
symmetry of S A (x;y) between u>2 and 5 (or p and q). 

B.2 Case of type BC 

We now consider the difference operator E^ S ' K ' of type BC m , defined by (|2.7j) . (|2.10j) . (I2.11|) 
with the function [x] of (|B,1|) . For the variables x = (x±, . . . , x m ) and the parameters \i = 
(f/,1, . . . , /is); k, we will use the multiplicative expressions together by setting 

z i = e(x i /uj 1 ) (i = l,...,m), u k = e(n k /u 1 ) (k = 1, . . . , 8), t = e(n/u 1 ). 
We first remark that, by (|2.13p . the difference operator can be rewritten in the form 

3 



f)-'- 

Ki<m: e=± r=0 
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Here the coefficients Af(x] fj,\8, k) (i = 1, ...,m), A®(x; fi\5, k) (r = 0,1,2,3) are obtained as 
follows by modifying the corresponding coefficients of Ex K 



A + (x- u\6 k) = ^ 17 K-[Xi ± Xj + K--r) 



l<j<m;j^i v 1 J 

A~(x;fi\5, k) = Af(-x; fj,\8, k) (i = 1, . . . , to), 



8 



J! e(u k q h r ;p) m , I ±i n 
r *Sl TT I j ' fr-Q19^ 

c " WipMPrhp) f = \ e( q -h r zf; P ) { ~ ' ' ' } 



with constants 

c r = e(u r /2uji), C r 



y-%T~ Vr (w r - 28 + mn + \ ^ W j (r = 0, 1, 2, 3). 



In the multiplicative expression, these constants are expressed as 

i _ i 

C = 1, Cl = -1, C 2 =p2, C 3 = -p 2; 

Co = d = i, c 2 = t m p-V 2 («i---«8)3, c 3 = t- m pV(«i---«8)-5. 

In view of formula (|B.2j) . we take a nonzero meromorphic function -P(x) satisfying the system 
of difference equations 

Tl(P(x))=qzfP(x) (i = l,...,m). 

m 

The simplest choice for such a function is given by P(x) = e(^j £ xfj. Then formula IpOl) 
implies that the difference operator E X ^ S ' K ' is expressed in the form 

e (p\S,k) = t - m+ i pq ^ ui . ..ugj-ip^gMWp^)-^ (B.3) 

as a constant multiple of the conjugation of a difference operator 

3 

l<i<m; e=± r=0 

by P(s). 

We also remark that the operator £^ S ' K ' has the following property with respect to shifting k 
to k — u>2 (t to t/p): 

£046,*) = ^qt-^^p^^gMS^-^p^ym+l^ t = e ( K / Wl ). 
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Together with (|B.3p . this also implies 

E (H\S, K ) = t -3(m-l) p 2m-l q m+^ Ui . . . Ug )-± p( x )™ £ (n\S,K-u, 2 ) p( x ym _ 

Setting k = X + L02 (t = sp, s = e(\/u>i)), we rewrite this formula as 

E (p,\S,X-+uJ2) = s -3(m-l) p -m+2 q m+i( Ul . . . U8 ykp( X ) m S^ S '^P(x)- m . (B.4) 

The difference operator Sx ^ appearing in (jB.4[) is essentially the same as the operator 
ll A + (h,n;x)" of type BC m defined by Ruijsenaars [23 I, (4.1)-(4.3)] (see also (3.1)-(3.9)); the 
difference between the two is only by an additive constant. In fact, we have 

3 

«A+(h,ti;xy =4^' A >-^t?(£ r ,...,£ r ;/45,A), (B.5) 

r=0 

where £0 — £2 = wi/2, £1 = £3 = 0, under the identification of the parameters 

"h k » =ii k -±(u 2 + S) (A: = 0,1,..., 7), "/i" = A, 
with fiQ = fig. 

By Theorem 12.31 we know that the kernel function of (|2.15[) 

m n 

$Bc(x;y\5,K) = mi Yl G(e 1 x j + e 2 yi + ^(5-K)\5) 

j=U=l ei,C2=± 

in the variables x = (x±, . . . , x m ) and y = (yi, . . . , j/„) satisfies the difference equation 

E^$ BC (x-y\8,K) = E^<S> BC (x;y\5,K), 
Jl k = ^(5 + k) - fi k (fc = l,...,8) 

under the balancing condition 

8 

{m-n-l)K-5+^^2fi k = (t mr - n - 1 q- 1 (u 1 ---u s )* = 1). (B.6) 
fe=i 

With G(x\S) = G+(x\5), (12. 4j) . this kernel function is given explicitly by 
$Bc(x;y|<5,«)=e(^)nn II r^sH^uf 

J=li=l ei,E2=± 

m n 

where wi = e(yi/uj\) (1 = 1,..., n), and f(x; y) = n ^ x 2 + y 2 + ^(k 2 — 5 2 ). We express 

j=i 1=1 

^Bc{x',y\S, k) in the form 

3> BC (x; y\5, k) = const • P(x) n P(y) m S BC (x; y), 

where 
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s B c{x-,y) = nn n ^^p,?). 

3=11=1 ei,e2=± 

It should be noted that, by the substitution k = A + u>2 (t = sp, s = e(A/o»i)), the function 

m n m n 

s BC (x;y) = Hll J] r( P !?5*"^< J ;j',?) = niI II Qfaj + em-l*) 

j=li=l ei,€2=± j=li=l ei,e2=± 

becomes manifestly invariant under the action of the Weyl group of type BC m ( resp. BCn) on 
the variables x = (x±, . . . , x m ) (resp. y = (yi, ... , y n )), and symmetric with respect to and (5 
(p and g). This function S B c{x; v) can be- thought of as the (BC m , -BC n )-version of Ruijsenaars' 
kernel function u S(x;y)" in [23J. 

We now have the functional equation 

E^P(x) n P(y) m S BC (x;y) = Ef\ 5 ^P{x) n P{y) m S BC {x;y) 

under the balancing condition (IB.6j) . Hence, by applying (|B.4p with k = A + 0J2, we obtain 

const • P{x) m 8^^P{x) n - m P{y) m S BC {x]y) 

= const • P{y) n £^P{x) n P{y) m - n S BC {x-y), 

namely, 

const • £^P(x) n - m P(y) m - n S BC (x;y) 

= const • £^P(x) n - m P(yr- n S B c(x;y). 

The constants in front of the both sides simplifies by (|B.6|) . to imply the functional equation 

(pq/s) m £^ s ^P(x) n - m P(y) m - n S BC (x;y) 

= {pq/s) n £^P{xT- m P{yr- n S BC {x; y), 
under the balancing condition 

8 

U2lA k = 8-(m-n-l)(\ + u 2 ) ((ui • • • u 8 )5 = q(p S r m+n+1 ). 

k=l 

In particular, when m = n we have 

£^S BC (x; y) = £^S BC (x; y), (B.7) 
under the balancing condition 

8 

I^Mifc = A + w 2 + <5 ((ui • • -us) 5 = spg). (B.8) 
fc=i 

The functional equation (|B.7|) is precisely the formula (4.24) of [231 I, Proposition 4.1], for 
"<5 = +" with the parameter 

"h k " =fj, k -±(u 2 + 5) (A; = 0,1,-.., 7), >" = A. 

7 

The balancing condition ()B.8|) corresponds to condition (4.26), "/x = 2za + | ^ h k " . Also the 

fc=o 

parameters jl k = |(A + + 5) — /Ufc = 1, ... ,8) for the y variables are consistent with 
u {-jRh) k = 1~ia- h k " (jfe = 0, 1, . . . , 7). (The constant V+(/t)" on the right side of [Ml I, 
(4.24)], arises as the difference of constant terms of two operators, as indicated in (|B.5|) .) 
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C E r (z; a\t) and Hi(z; a\q, t) as interpolation polynomials 

In Section we presented some explicit expansion formulas for Koornwinder polynomials at- 
tached to single columns and to single rows, in terms of invariant Laurent polynomials E r {z; a\t) 
and Hi(z;a\q,t), respectively. These Laurent polynomials E r (z;a\t) and Hi(z;a\q,t) are essen- 
tially the same objects as the SC m -interpolation polynomials of Okounkov [18] attached to 
single columns and single rows. Our Theorems 15. II and 15.21 provides with explicit expressions for 
the corresponding special cases of the binomial expansion of Koornwinder polynomials in terms 
of £>C m -interpolation polynomials as discussed in Okounkov [18] and Rains [19] . 

In the notation of Section [5l the BC m interpolation polynomial P£(z;q,t,a) attached a par- 
tition A = (Ai, . . . , A m ) [18] is characterized uniquely, up to multiplication by a constant, as a 
Laurent polynomial in HC[z =l: ] of degree |A| satisfying the following conditions: 

(0) P£(z;q,t,a) is ^-invariant in the shifted variables zt s a = {z\t m ~ l a, Z2t m ~ 2 a, . . . , z m a), 
where 5 = (m — 1, m — 2 , 0). 

(1) P^(q^] q, t, a) = for any partition \i = . . . , fj, m ) such that fj, A. 



In this section we show that the Laurent polynomials E r (zt a;a\t) and Hi(zt a;a\q,t) co- 
incide, up to constant multiplication, with the interpolation polynomials P,* ^ (z\ q, t, a) and 
Pj!jJz;q, t, a), respectively. We prove that these polynomials actually have the interpolation 
properties as mentioned above. 

Proposition C.l. For each r = 0, 1, . . . , m, let E r (z; a\t) be the W -invariant Laurent polyno- 
mial defined in (|5.ip . 



Proof. Since E r {z; a\t) is symmetric in z = (z±, . . . , z m ), from the recurrence relation (|5.5|) we 



(2) P A % 



q,t,a) / 0. 




have 



E r (z;a\t) = E r (zi, . . . , z m -i;ta\t) + E r -\{z\, . . . , z m -i; a\t)(z m ; a) 



and hence by setting z m = a, 



E r (zi 



z m -i,a;a\t) = E r (z 1 , . 



z m -i;ta\t). 



This implies 



E r (zi 



.m—r—l 



a, . . . ,a;a 






z r ;t m - r a\t) 



for r = 0, 1, . 



, m, which proves (0). In particular, we have 



E r (z u ..., z r -i,t m - r a, t™-*- 1 ^ ...,a; a\t) = 0, 

which implies E r (q fJ, t s a;a\t) = for any partition such that l(fi) < r, namely, \i ~£> (l r ). State- 
ment (0) also implies 



E r (gr^a, . . . , qt m ~ r a, t m " r - 1 a, . . . , a; a\t) = (-If (t m ~ r a; qt m - r a)t, r - ■ 
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From this proposition we see that E r (zt a;a\t) is a constant multiple of the interpolation 
polynomial P? lr Jz; q, t, a) for each r = 0, 1, . . . , m. 

Proposition C.2. For I = 0,1,2,..., let Hi(z;a\q,t) be the W -invariant Laurent polynomial 
defined in (|5.3p . 

(1) For any partition \i = (jii, . . . , fi m ) with \l\ < I, i.e., /u 7$ (I), Hi(q^t s a; a\q, t) = 0. 

(2) H l (g l t m - 1 a,t m - 2 a ) ...,a;a\q,t) = (t) g ,,(g , t 2 ( m - 1 )o 3 ),,,. 

Proof. By the definition (|5.3|) . the Laurent polynomials Hi(z;a\q,t) (I = 0, 1,2. . .) satisfy the 
recurrence formula 

1 (t) l-r 

Hi(z;a\q,t) = y~] - ^ r H r (zi, . . . , z m -\\ a\q,t)(z m ; q r t m ~ 1 a} g: i- r . 

r=0 \Q)q,l-r 

Since these Laurent polynomials are symmetric in z = (z±, . . . , z m ) (Lemma I5.4p . we also have 

(t) i 

Hi(z; a\q, t) = ^ q ' " r (zr, t m ~ 1 q r a) q ^ r H r (z 2 , ■ ■ ■ , z m ; a\q, t). 
r= Q \*/ qJ~ r 

We prove proposition by the induction on m. For an arbitrary partition = (fii, . . . ,/i m ), we 
specialize this recurrence formula to z = q^t s a: 

miqH'a; a\q, t) = ^j^^H^a- qH^^^H^H^a, q^a; a\q, t). (C.l) 

By the induction hypothesis we have 

H r (q^t m ~ 2 a, . . . ,q» m a;a\q,t) =0 
for r > n\. For r < fj,±, we see 

{q^t m ~ l a- q r t m ~ l a) q ^ r = (-l) l ~ r (q-^+ r ) q>l _ r (q^+ r t 2 ^a 2 ) q ^ r . 

Since — /xi + r < 0, (q~ ,J ' 1+r )q ) i-r = if / — r > fj,\ — r, i.e., \i\ < I. This means that all the terms 
on the right side of (jC.ip vanish when [i\ < I, and hence, Hi{q^t & a\a\q,t) = 0. When \i = (/), 
the only nontrivial term arises from r = 0: 

H l {q l f n - 1 a,t m - 2 a,...,a-a\q,t) = (qk™' 1 a; t m ~ l a) = (t^qH 2 ^- 1 ^ 2 )^. U 

This proposition implies that Hi(zt s a;a\a,t) is a constant multiple of the interpolation poly- 
nomial P(n(z\ q, t, a) for each I = 0, 1, 2, 
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